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Abstract. Let r ; 5 x 5 — > IR+ be the jump rates of an irreducible random 
walk on a finite set S, reversible with respect to some probability measure m. 
For a > 1, let g : N ^ R+ be given by g(0) = 0, g(l) = 1, g(k) = (k/k - 1)", 
fc > 2. Consider a zero range process on S in which a particle jumps from a 
site X, occupied by k particles, to a site y at rate g{k)r{x, y). Let stand for 
the total number of particles. In the stationary state, as Af t c», all particles 
but a finite number accumulate on one single site. We show in this article 
that in the time scale N^^°' the site which concentrates almost all particles 
evolves as a random walk on S whose transition rates are proportional to the 
capacities of the underlying random walk. 



1. Introduction 

Fix a finite state space S of cardinality k > 2 and consider an irreducible con- 
tinuous time random walk {Xt : t > 0} on 5* which jumps from a; to y at some rate 
r{x,y). Assume that this dynamics is reversible with respect to some probability 
measure m on S: m{x)r{x, y) = m{y)r{y, x), x, y € S. Denote by cap^ the capacity 
associated to this random walk: For two disjoint proper subsets A, B of S, 

caps(A,B) = inf i ^ m(x) r{x,y) {f{y) ~ f{x)}\ 

where B{A, B) stands for the set of functions / : 5 ^ K equal to 1 at A and equal 
to at B. When A ~ {a;}, B — {y}, we represent cap5(A, B) by cap5(a;,?/). 

Let Mi, be the maximum value of the probability measure m: Mi, = max{TO(a;) : 
x G S'} and denote by Si, the sites where m attains its maximum value: Si, = 
{x £ S : m{x) = Mi,}. Of course, in the symmetric, nearest-neighbor case, where 
r{x, y) = 1 if y = X ± 1, modulo k, and r{x, y) = otherwise, m is constant and Si, 
and S coincide. 

Fix a real number a > I. Let g : N ^ M be given by 

(n \ ^ 
' "-2' 

so that nr=i = IT-"', n > 1. Consider the zero range process on S in which a 
particle jumps from a site x, occupied by k particles, to a site y at rate g{k)r{x, y). 
Since g is decreasing, the dynamics is attractive in the sense that particles on sites 
with a large number of particles leave them at a slower rate than particles on sites 
with a small number of particles. 

The total number of particles is conserved by the dynamics, and for each fixed 
integer > 1 the process restricted to the set of configurations with N particles. 
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denoted by Em, is irreducible. Let he the unique invariant probability measure 
on En- When a > 2, the measure exhibits a very peculiar structure called 
condensation in the physics literature. Mathematically, this means that under the 
stationary state, above a certain critical density, as the total number of particles 

I oo, only a finite number of particles are located on the sites which do not 
contain the largest number of particles. 

Condensation has been observed and investigated in shaken granular systems, 
growing and rewiring networks, traffic flows and wealth condensation in macroeco- 
nomics. We refer to the recent review by Evans and Hanney [4J. 

Several aspects of the condensation phenomenon for zero range dynamics have 
been examined. Let the condensate be the site with the maximal occupancy. Precise 
estimates on the number of particles at the condensated, as well as its fluctuations, 
have been obtained in [HI [H |5j. The equivalence of ensembles has been proved 
by GroBkinsky, Schiitz and Spohn [5]. Ferrari, Landim and Sisko [B] proved that 
if the number of sites is kept fixed, as the total number of particles N ^ oo, the 
distribution of particles outside the condensated converges to the grand canonical 
distribution with critical density. Armendariz and Loulakis [I] generalized this 
result showing that if the number of sites k grows with the number of particles N 
in such a way that the density N/k converges to a value greater than the critical 
density, the distribution of the particles outside the condensate converges to the 
grand canonical distribution with critical density. 

We investigate in this article the dynamical aspects of the condensation phe- 
nomenon. Fix an initial configuration with the majority of particles located at one 
site. Denote by X^^ the position of the condensated at time i > 0. In case of ties, 
X^^ remains in the last position. The process {X^ ■ t > 0} evolves randomly on S 
according to some non-Markovian dynamics. 

The main result of this article states that, for a > 1, in the time scale iV^+", 
the process {X^ ■ t > 0} evolves asymptotically according to a random walk on 5*^, 
which jumps from x to y at a rate proportional to the capacity cap_5(a;, y). In the 
terminology of |2j| , we are proving that the condensate exhibits a tunneling behavior 
in the time scale A^^+". 

This article leaves two interesting open questions. The techniques used here rely 
strongly on the reversibility of the process. It is quite natural to examine the same 
problem for asymmetric zero range processes where new techniques are required. 
On the other hand, the number of sites is kept fixed. It is quite tempting to let 
the number of sites grow with the number of particles. In this case, in the nearest 
neighbor, symmetric model, for instance, the condensate jumps from one site to 
another at rate proportional to the inverse of the distance. The rates are therefore 
not summable and it is not clear if a scaling limit exists. 

Simulations for the evolution of the condensated have been performed by Go- 
dreche and Luck [7]. The authors predicted the time scale, obtained here, in which 
the condensate evolves and claimed that the time scale should be the same for non 
reversible dynamics. 



2. Notation and results 

Throughout this article we fix a finite set S of cardinality k > 2 and a real 
number a > 1. For each So 'Z S consider the set of configurations En^Sq: N > 1, 
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given by 

En,So ■■= {V e : ^ 7?, = iV} , 

xeSa 

where N = {0, 1, 2, ...}. When 5*0 = S, we use the shorthand for En,s- Define 
a{n) = for n > 1 and set o(0) = 1. Let us also define g -.N ^ , 

5(0) = 0, ff(l) = l and g(n) = , ^ \. , n>2, 

a[n — 1) 

in such a way that nr=i ~ '^C'^)) n > 1, and {g{n) : n > 2} is a strictly 
decreasing sequence converging to 1 as n t oo. 

Consider a random walk on S with jump rates denoted by r. Its generator Cs 
acts on functions / : 5* ^ M as 

{Csf){x) = Y,r{x,y){f{y)-f{x)}. 
yes 

Assume that this Markov process is irreducible and reversible with respect to some 
probability measure m on S: 

m(x)r(x,y) = m{y)r{y,x) , x,y € S . (2-1) 

Let be the maximum value of the probability measure m, let 5* C S be the 
sites of S where m attains its maximum and let be the cardinality of 5*: 

= max{m(a;) : x G S} , = {x G S : m{x) = M*} and = l^*! . 

In addition, let m*(a;) = m{x)/Mi, so that m*(a;) = 1 for any x G /S*. Denote by 
Ds the Dirichlet form associated to the random walk: 

Dsif) = ^ E m{x)r{x,y){f{y)-f{x)y (2.2) 

x,yGS 

for / : 5 — + M, and denote by cap5(a;,y) the capacity between two different points 
X y G S' 

caps(a;,y) = inf {Ds{f)} , (2.3) 

feB{x,y) 

where the infimum is carried over the set B{x, y) of all functions f : S ^ M. such 
that f{x) = 1 and /(y) = 0. 

For each pair x,y G S,x ^ y, and rj G En such that % > 0, denote by a^^r/ the 
configuration obtained from r/ by moving a particle from x to y: 

{rjx — I ioi z = X 
r)y + l toi z = y 
riz otherwise . 

For each N > 1, consider the zero range process defined as the Markov process 
{ri^{t) : t > 0} on E^ whose generator Ljv acts on functions F : Ej^ — > M as 

{L^F){r^) = M r{x, y) {i^(<7-^77) - F{r,)) . (2.4) 

x^y 

The Markov process corresponding to Ljv, > 1, is irreducible and reversible with 
respect to its unique invariant measure given by 

Atjv(r? = M — -— := — , tiGEn , 

Zn,s j--^ a(?7x) Zn,s a{r]) 
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where, for any So S, 

and Zn.So normalizing constant 

In Section [3] we show that the sequence {Z^^s '■ N > 1} converges as TV t oo. 
This explains the factor iV" in its definition. The precise statement is as follows. 
For a; in 5* and /t > 2, let 



so that T{a) = Tx for any x e 5*, and define 



Proposition 2.1. For every k> 2, 

hm = Zs ■ 

Denote by Dn the Dirichlct form associated to the generator L^. An elementary 
computation shows that 

^A^(^) = E f^N{v)9iVx)rix,y){F{a^yv)-F{r,)r , 

a:. yes TJGEn 

for every F : En K. 

For every two disjoint subsets A, B of Em denote by Cn{A, B) the set of functions 
F -.En ^M. defined by 

Cn{A, B) {F : F{rj) = 1 V G A and F{() = V ^ G 5} . 

The capacity corresponding to this pair of disjoint subsets A, B is defined as 

capjv(^,S) := mi{DN{F) : F eCNiA,B)} . 

Since Dn{F) =Dn{1~F), cap^(A,B) =cap^(B,A). 

Fix a sequence {£n : > 1} such that 1 <^ £n N and, for each z e S \ 5*, 
fix a sequence {bN{z) : ^ > 1} such that 1 ^ bN{z): 

lim ^jv = 00 , lim £n/N — , and lim bN{z) = 00 , (2-6) 
for all z e S \ Si,. For x in 5^, , let 

£-N - {v(^EN:Vx>N-iN,Vz< bN{z) , z ^ . 

Obviously, £^ 7^ for all a; G 5* and every N large enough. In the case where the 
measure m is uniform, the second condition is meaningless and the set £^ becomes 
t'i, = {v&EN ■.Vx>N -In}- 

Condition £n/N ^ is required to guarantee that on each set £^ the proportion 
of particles at cc G S**, i.e. rj^/N, is almost one. As a consequence, for N sufficiently 
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large, the subsets x € Si,, are pairwise disjoint. ^From now on, we assume that 
N is large enough so that the partition 

-£^iv = ( U £?v) U (2.7) 
is well defined, where A^r is the set of configurations which do not belong to any 

set ^ X ^ • 

The assumptions that t oo and that bj^{z) ] oo for all z ^ Si, arc sufficient 
to prove that ij,n{'^n) ^ 0, as we shall see in p.2p . and to deduce the limit of 
the capacities stated in Theorem 12.21 below . In particular, in these two statements 
we may set bN{z) — N, z ^ Si,, in order to discard the second restriction in the 
definition of the sets x E Si,. We need, however, further restrictions on the 
growth of and b]\[{z) to prove the tunneling behaviour of the zero range processes 
presented in Theorem 12.41 below. 

To state the first main result of this article, for any nonempty subset S'^ of Si,, 
let £^(5"^) = U^esi£-Ni and let 

1 

u°'{l-u)"du. (2.8) 

Theorem 2.2. Assume that Ki, > 2. Fix a nonempty subset Sl C Si, and denote 
S^ ^ Si,\ Si- Then, 

hm 7Vi+"cap^(£^(^i),£^(52)) ^ ^ 1^ ^ cap5(x,y). 

Note that the right hand side depends on the sites not in Si, through the capac- 
ities ca.pg{x,y). In the case where the measure m is constant, 5*^, = S, Mi, = k^^ 
and the right hand side becomes 

To prove Theorem [221 we derive a lower and an upper bound for the capacity. In 
the first part, we need to obtain a lower bound for the Dirichlet form of functions in 
Cn{£-n{SI), £Ar(5'^)). To our advantage, since it is a lower bound, we may neglect 
some bonds in the Dirichlet form we believe to be irrelevant. On the other hand, and 
this is the main difficulty, the estimate must be uniform over Cn{£.n{SI), £.n{SI))- 
As we shall see in Section^ the proof of a sharp lower bound gives a clear indication 
of the qualitative behavior of the function which solves the variational problem 
appearing in the definition of the capacity. With this information, we may propose 
a candidate for the upper bound. Here, in contrast with the first part, we have 
to estimate the Dirichlet form of a specific function, our elected candidate, but we 
need to estimate all the Dirichlet form and can not neglect any bond. 

For each ry g E^, let stand for the probability on the path space I?(M+, i?jv) 
induced by the zero range process {??^(i) '■ t > 0} introduced in (|2.4|) starting from 
T] G En. Expectation with respect to is denoted by E^. In addition, for any 
A C En, let Ta denote the hitting time of A: 

Ta ■■= inf > : tj^ (t) e A} . 
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Remark 2.3. It is well known (see e.g. Lemma 6.4 in ^) that the solution of the 
variational problem for the capacity is given by 



N 



The candidate proposed in the proof of the upper hound provides, therefore, an 
approximation, in the Dirichlet sense, of the function ¥31^2. 

The second main result of this article states that the zero range process exhibits 
a metastable behavior as defined in [2]. Fix a nonempty subset A of E^. For each 
t > 0, let T/^ be the time spent by the zero range process {ri^ {t) : i > 0} on the 
set A in the time interval [Q,t]: 

f\{v''{s)eA}ds 
Jo 

and let be the generalized inverse of T/^ : 

Sf := sup{s > : r/ < . 

It is well known that the process {ri^'"^{t) : t > 0} defined by ri^'^{t) = rj^ (Sf) is 
a strong Markov process with state space A [5] . This Markov process is called the 
trace of the Markov process {ri^{t) : t > 0} on A. 

Consider the trace of {r]^ (t) : i > 0} on Eat (5**), referred to as 77^" (t). Let 
^jv : ^NiSi,) I— > Si, be given by 

and let := *Ar(r;^« (i)). 

We prove in Theorem 12.41 below that the speeded up non-Markovian process 



{X^^^i : t > 0} converges to the random walk {Xt : t > 0} on 5* whose generator 
Lg^ is given by 

(Ls./)(x) = ] c^P^^(^' y) {/(y) - /(^)} • (2.9) 

Ah T [a) la ^ 

y to* 

For X in 5*, denote by Pj. the probability measure on the path space -D(M+,S'^) 
induced by the random walk {X^ : t > 0} starting from x. 

Theorem 2.4. Assume that > 2. // (|2.6[) holds and 

.1+q(k-1) 

lim n mJz)-'"'^'^ = (2.10) 

then, for each x G Si,, 
(Ml) We have 

(M2) For any sequence G N > 1, the law of the stochastic process 

: t>Q} under converges to as TV "f 00; 
(M3) For every T > 0, 

i-T 



N 



lim sup E, 



/ l{ry^(siV"+i) e Ajvjds 
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If K > K*, in order to fulfill conditions (|2.6|) and ()2.I0p . we can take, for instance, 

= ~ fovzeS\S, and £n = N^^^^-^^ 

log(m*(z)) 

if > 3, = ivi/[«-(i/2)l if = 2. 

According to the terminology introduced in [5] , Theorem 12.41 states that the 
sequence of zero range processes {rj^ {t) ■ t > 0} exhibits a tunneling behaviour on 
the time-scale iV"+^ with metastates given by {£^ : a; e S**} and limit given by 
the random walk {Xt : t > 0}. 

Property (M3) states that, outside a time set of order smaller than N"^^, one 
of the sites in 5* is occupied by at least N ~ £n particles. Property (M2) describes 
the time-evolution on the scale N°'^^ of the site concentrating the largest number 
of particles. It evolves asymptotically as a Markov process on 5* which jumps from 
a site X to y at a rate proportional to the capacity ceipg{x,y) of the underlying 
random walk. Property (Ml) guarantees that the process starting in a metastate 
£^ thermalizes therein before reaching any other metastate. 

Remark 2.5. In [Q, it is shown that, in the case the number of sites increases with 
the number of particles, the highest occupied site contains a nonzero fraction of the 
particles in the system. This result includes the case 1 < a < 2. In contrast, when 
the number of sites is kept fixed, it seems to have been unnoticed in the literature 
that the condensation phenomenon appears also for 1 < a < 2. More precisely, if 
K ^jv < then 

lim fJ-NiVx > N - In) = , Vx e S** . 

Moreover, given that particles concentrate on x ^ S^,, the distribution of the configu- 
ration on S\{x} is asymptotically given by the grand- canonical measure determined 
by rrii,: For any x in Si,, 

1 m^(z)^- 



lim sup 



1^n{vz = Cz , z x\t]x > N - £n) - Y\_ 



= 



where := {C € N■5\^^l : („, < tj^}. There is just a small difference between 
the cases 1 < a < 2 and a > 2. While in the former, the variables {riz ■ z € S^} 
do not have finite expectation under the critical grand- canonical measure, they do 
have finite expectation in the latter case. 

In 13], we have proved Theorems 12.21 and 12.41 in the case where the rates r(-, •) 
in the definition of the sequence of zero range processes corresponds to a random 
walk on a finite complete graph. Since it covers the case k = 2, we may suppose 
that K > 3. 

3. The stationary measure 

In this section, we prove Proposition l2.1l The proof relies on four lemmata. We 
first show that the sequence Zj^^s is bounded below by a strictly positive constant 
and above by a finite constant. Let 

1 

'■'N,i 



•^-^ a(r]) 



and note that Zn^s < ^Ar,|S|- 
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Lemma 3.1. For each k > 2, there exists a constant > 0, which only depends 
on a and k, such that 

1 < Zn,S < -^iV.re < Ak . 

Proof. Choose x in S-^, and denote by ^ the configuration in Ej\j such that ^(x) — N, 
= for y 7^ x. By definition, m*(x) = 1 so that ^jv,s > N°'mi/a{^) = 1, 
which proves the lower bound. 

We proceed by induction to prove the upper bound. The estimate clearly holds 
for K = 2. Assume that it is in force for 2 < k < A:. The identity 

1 Z 
= ^1]v^ + E aijHN-j)} 

permits to extend it to «; = fc. □ 

For any £ > 1, let En,s{^) be the subset of -Ejv,s of all configurations with at 
most N — i particles per site: 

EN,s{e) = {?? e En,s :Vx<N-e,yxeS} . 

Next lemma shows that the measure /it at is concentrated on configurations in which 
all particles but a finite number accumulate at one site. 

Lemma 3.2. There exists a constant > which only depends on a and k, such 
that for every integer £ > 0, 

Proof. We proceed by induction on k. For k = 2 the statement is easily checked. 
Now, suppose the claim holds for 2 < k < k — 1. Fix some x in S. The left hand 
side of the inequality in the statement can be written as 



This sum is equal to 



o<i<e/2 e/2<i<N~e ^ ' ^ ' ce-EN-,,s\{x}(^-i) 

where the second sum is equal to zero if {i : £/2 < i < A'^ — is empty. We examine 
the two terms of this expression separately. By the induction assumption, the first 
sum is bounded above by 

t A— 1 

a{i)a{N-i) {i-i)"-^ ' 
By the previous lemma, this sum is less than or equal to 

2"-iCfc-i ^ A" ^ 2"-i(7fc-iZjv,2 ^ 2«-iCfe_iA2 



.0 ^HHN-i) 



lla—l — (<x—\ 
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On the other hand, by Lemma 13.11 and the induction assumption for k = 2, the 
second term in is less than or equal to 



^-^ a{i)a(N — i) 
e/2<i<N~e ^ ' ^ > 

This concludes the proof of the lemma. □ 

For > 2, < £ < A^, X € S", denote by i?^'^ the set of configurations in E-^,s 
with at least N ~ I particles at site x: 

Recall the definition of the set Si,. Next lemma shows that the /i at -measure of the 



set E^jf decays exponentially if x does not belong to Si,. 



Lemma 3.3. For each k>2, there exists a finite constant C„, depending only on 
K and a, such that 

for all N > e. 

Proof. Fix K > 2 and x in 5*. The expression on the left hand side of the statement 
of the lemma is bounded by 

By Lemma |3. II and since mi,{x) < 1, the second term is less than or equal to 



N-l 



which concludes the proof of the lemma. □ 
If £ < N/2, the sets {£;^' : x e S} are pairwise disjoint and 
En,s\ U E^/ = EN^s{i+l) . 

It follows from the two previous lemmata that the sum in the definition of Zpf^s 
restricted to configurations on E^^, x £ Si,, is close to k^^Zn.s for £ (and conse- 
quently N) large. 

Lemma 3.4. For each k > 2, there exists a constant > 0, which only depends 
on a and n, such that for every integer ^ > and x £ Si,, 



sup 

N>21 



ml Zn,s 



< 



Proof. As we have observed, for < £ < N/2, 



71 7] 
, Txrrv ^* 



a(ri) ^ air/) 
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By symmetry, for x, y in 5*, 

\ - mj _ X - m^; 



Hence, if x belongs to S'^,, 



TJ 77 77 

ml , "i; ^ \ - 



The statement now follows from the two previous lemmata. □ 
We are now in a position to prove the main result of this section. 



Proof of Provosition \2.1\ Fix a site a; in S**. By the previous lemma, 
lim h-^Zns = lim lim N" V 

Since x belongs to S*^, the previous sum is equal to 

As iV I cx) and i ^ oo, this expression converges to 

^ \^ JiA. - TJSrUhML - TTr 
^ ^ a(0 " a(j) " -Li'^- 

This concludes the proof of the proposition. □ 

We close this section showing that 

lim A*iv(AAf) = . (3.2) 

N^oo 

Recall the definition of the set 5'^ and of the sets x E S^. Since 

Ajv = [ U {v:Vz>bN{z)}]\J [f] {v:v.<N-eN} 
z£S\St xes^ 

intersecting the second set with the partition A = r]zes\St, {jl '■ < N ~ ^jv} and 
A'^, we get that 

Aat C U E'/- y En,s{^n + 1) , 
zes\s^ 

where cat — min{^jv,7V— hN{z) : z £ S \ Si,}. Hence, assertion (13.21) follows from 
Lemma 13.21 assumption (|2.6p and Lemma [ 
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4. Lower bound 

In this section we prove a lower bound for the capacity. It might be simpler in 
a first reading to assume that m is constant so that S — S^,. 

For b,£ > 3 and x, y in Si,, x y, consider the tube defined by 

L^n"" - {v(^EN:Vx+Vy>N~£;r],<b,zeS\S.y 

Clearly, = for any x,y £ 5*. We claim that for each x £ Si, and every TV 
sufficiently large 

L^^nL^" C £^ , y,zeSA{^}- (4.1) 

Indeed, let ry £ n . First, b < int z£s\St bniz) for any N sufficiently large 
in view of (|2.6p . On the other hand, rjz < i because rj belongs to . Hence, 
Vx > N — 2£ since rj belongs to L^^. Since ijy ^ oo, this shows that r]x > N — i^, 
for N large enough and we conclude that rj £ Moreover, it follows from this 
argument that, for N sufficiently large, 

n i^™ 7^ if and only if {x, y} n{z,w} ^0 . (4.2) 

Proposition 4.1. Assume that > 2. Fix a nonempty subset S^ C S^, and denote 
S'^ = Si,\Sl Then, 

lnninfiVi+"cap^(£^(5,i),£A,(52)) > \—— ^ caps(x,2/). 

Proof. Fix a function F in CAr(£jv(S'^), £.n{S^)). By definition, 

^^(^) = ^ E E MA'(r/)r(z,u;)g(77.){F(a^"'ry)-F(r;)}2. 

We may bound from below the Dirichlct form Dn{F) by 

^E E E E Mv)riz,w)girj.){Fia^--rj)~^Firj)r . 

In this inequality, we are neglecting several terms corresponding to configurations 
•q which do not belong to ^x^s'^,y£S'^L^N' . On the other hand, some configurations 
are counted more than once because the sets {L^^^ : x £ Sl,y £ S^} are not 
disjoints. However, by (|4.2p . if and L^'^ are different strips and rj belongs 

to n i^'^ then x — x' and y ^ y' (recall that = Lr'ff). In consequence, 
Vx > N — 2£. In particular, for TV large enough, ij and cr^™?? belong to for all 
z,w £ S, so that F{a^'^r]) — F{rf) because F is constant on £^. 

The proof of the lower bound has two steps. We first use the underlying random 
walk to estimate the Dirichlet form Dn{F) by the capacity of this random walk 
multiplied by the Dirichlet form of a zero range process on two sites. This remaining 
Dirichlet form is easily bounded by explicit computations. 

Fix X £ Sl, y £ S^. Denote by O^, x £ S, the configuration with one and 
only one particle at x, and agree that summation of configurations is performed 
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componentwise. The change of variables £, = rj — dz shows that 



1 ]\ia ™4 

9 E E 7 j^rn4z)r{z,w){F{^ + ^^)-F{( + ^z)V 



2 .t^s ^^.^ 
This sum is clearly bounded below by 



1 A/'" rr? * 

9 E E y -^m4z)r{z,w){F{^ + :,^)~F{^ + Vz)r ■ 

5^<6-1,V2GS\S. 

Fix a configuration ^ in E^-i and consider the function / : S* ^ R given by 
/(t;) = + - + + 0.) - F(e + 0,;)}. Note that f{x) = 1, 

/(j/) = 0. Moreover, if we recall the expression \2.2\ of the Dirichlet form of the 
underlying random walk, 

i ^ m,{z)r{z,w){F{^ + l}^)~F{i + liz)Y 

z,w£S 

= ■^^Dsif){FiC + ^.)-Fi^ + Vy)r . 
Since f{x) — 1, f{y) = 0, the previous expression is bounded below by 



-V caps(x, y) {F{^ + 0,) - F{^ + dy)} 
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Up to this point we proved that the Dirichlet form of F is bounded below by 

^j<b-i,Vzes\s. 

Fix xo e si, yo G 52 and let := 5" \ {2:0, Vo}- For each > 0, let Bk = 3^°'^° 
be the set of configurations on 6*0 given by 

Bk = { C e : ^ C« = fc ; G < 6 ~ 1 , ^ e -5 \ 5. } . 

veSo 

For C in Bfc, let G(; : {0, . . . , iV - 1 - fc} ^ R be defined as Gi;{i) = F{^), where ^ e 
i?jv- 1 is the configuration given hy — Cv, v & Sq, S,xo = * and £,yg — N~l~k~i. 
With this notation, for xq E Sl, yo ^ S'^ fixed, we may rewrite the second sum in 
the previous formula as 

N \ ^ X ^ ml \ ^ 1 



Zns^^ a(C) ^ a(i)a(N -1-k-i) 

because TO*(a;o) = m^(?/o) = 1- Note that G,^ is equal to on the set {0, . . . , ^at — fc}, 
and equal to 1 on the set {N — £n, . . . ,N — I — k}. We may therefore restrict the 
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sum over i to a subset. It is easy to derive a lower bound for 

The function G which minimizes this expression is given by G{N ^ £n) — 1, 

G{i + 1) - G{i) = a(i) a{N - I - k - i) , i e [In - k , N - £n - l] , 

where Kj^i is a normahzing constant to ensure the boundary condition G{£n — k) = 
0. The respective lower bound is 

ENixo,yo) { E a{i) a{N - I - k - i)^ 

i—lN 

This expression depends on the configuration C only through its number of particles. 
Moreover, for every fixed fc, A^^+-^"SAr(xo, ?/o) converges to I^^ as N ] oo. 
In conclusion. 



By Proposition 12.11 and the above conclusions, as iV f oo, the right hand side 
converges to 

1 X ^ X ^ X ^ ml 



E caps(x,y) Y 



Recall that £ and 6 are free parameters introduced in the definition of the strip 
. Thus, letting b '\ oo and then £ ^ oo, the second sum in the last expression 
converges to 



V V — = TT V m^zy ^ TT r = 

^ ^ a(C) ^ a(j) J-i ^ K^r(a) 

fe>0 C6-Efc,s\{x,„} zCiS\{x,y} J>0 ' zeS\{x.v} * ^ ' 



For the last equation we have used the explicit formula of Zg presented just before 
Proposition This proves the lemma. □ 

5. Upper bound 

We prove in this section an upper bound for the capacity. As in the previous 
section, it might be simpler in a first reading to assume that m is constant so that 
S = S-k, 

Proposition 5.1. Assume that > 2. Fix a nonempty subset Sl C Si, and denote 
= Si,\ Si- Then, 

limsup7Vi+"cap^(£Ar(S']),£Ar(S'2)) < — ^ V caps(x,y). 

In view of the variational formula for the capacity, to obtain an upper bound 
for cap^(£Ar(5^), £jv('S'^)), we need to choose a suitable function belonging to 
^^(£^(5^), £Ar(5'^)) and to compute its Dirichlet form. Recalling the proof of 
the lower bound, we expect this candidate to depend on the function which solves 
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the variational problem for the capacity of the underlying random walk and on the 
optimal function for the zero range process with two sites. 

To introduce the candidate, iix x G S^, y G and recall the definition of 
the tube . In view of the proof of the lower bound, the optimal function 
F £ CNi£-NiSl),EN{S'i)) on the tube should satisfy 

+ + = {{^y{w)-{^y{z)}{Fi^ + d,)-F{^+<)y)} 

= {fxyH - fxy(^)} {G(e. + 1) - G(e.)} , 

where fx,y is the function which solves the variational problem (|2.3p in B{x, y) for 
the capacity of the underlying random walk, and G is the function appearing in the 
proof of the lower bound. 

Since, on the tube i^^, X^z/a y ^ and G is a smooth function, paying 
a small cost we may replace £,x in the previous formula by + X^zgA^z f*-"^ ^'^y 
suitable set ^ C S" \ {a:, j/}. The natural candidate on the strip L^j^^ is therefore 

K-l 

F.y{(,) ■■= Y.^i.y{z,) - fxy(^, + l)} G(6, + ' ' ' + C.,) , 

where x = zi, Z2, . . . , = y is an enumeration of S such that fxy(2:j) > fx,y(zj+i) 
for 1 < j < K. A simple computation shows that this function has the required 
properties listed in the previous paragraph. 

Since the tubes L^^, x e Si, y € Si, are essentially disjoints, the candidate F 
should be equal to Fxy on each tube L^^f and equal to some appropriate convex 
combination of these functions on the complement. 

We hope that this informal explanation helps to understand the rigorous and 
detailed definition of the candidate we now present. Let 2) C K'^ be the compact 
subset 

T) := {-It e K;^ : ^ = 1} ■ 

For each different sites x,y G S and (5 > 0, consider the subsets of D 

:= {u e Ti :ux > I - S} and iL^^ {u e Ti : Ux + Uy > 1 - S} 

Clearly L^^ = for any x,y G S. 

Fix an arbitrary < e < 1/6 and a; in 5. Let %^ = OC'f/e) := \ Df,, y =^ x. 
Since Xy, y G S \ {x}, is a collection of pairwise disjoint compact subsets of T>, 
there is a family of smooth functions 

e^:2)^[0,l], yGS\{x}, 

such that J2yeS\{x} Qy(^) — 1 ^o'' " 6y(u) — 1 for all u in %y and 

yGS\{x}. 

Clearly, the sets are macroscopic versions of the strips . The functions 
&y will be used to define the candidate function in the complement of the cylinders 

TX::y 

■ 

Let H : [0, 1] [0, 1] be the smooth function given by 

1 r<i>{t) 

H{t) := — / , 

la Jo 

where la is the constant defined in (|2.8p and (j) : [0, 1] [0, 1] is a smooth bijective 
function such that ^(t) + 0(1 - = 1 for every t G [0, 1] and (j){s) = Vs G [0, 3e]. 
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It can be easily checked that 

H{t)+H{l-t) = l , Vte[0,l], (5.1) 

7J|[Q 3^] = and -ff|[i_3£ ij = 1. The function iJ is a smooth approximation of the 
function G which appeared in the proof of the lower bound. 

Recall that a; g 5 is fixed. For each y G S \ {x} consider the function fxy : S 
[0, 1] in B{x, y) such that 

£'s(fxy) = cap5(a;,y) = inf Ds{f) ■ 

It is well known that fxy{z) is equal to the probability that the random walk with 
generator £5 reaches x before y when it starts from z. 
For each y (z S \ {x} fix an enumeration 

X ^ zi, Z2 , ... , = y (5.2) 

of S satisfying f^^iyi^j) > ixy{zj+i) for 1 < j < k — 1 and define Fr^y : Ejy ^ M as 
the convex linear combination 

K-l 

Pxyiv) ■■= ^{f^y{Z]) - fxy(^j + l)} Fly{l]) , rjeEN , 
J = l 

where each F^y : En ^ R, 1 < j < t — 1, is given by Fly{rj) = H{rix/N) and 

^ivW := ^f(^+niin{^^77.,;e}) , V & En , (5.3) 

1=2 

for 2 < j < K- 1. 

The function F^y just defined is a smooth approximation of the function F^.y 
defined at the beginning of this section. It is therefore the candidate to solve the 
variational problem for the capacity on the tube . It remains to define F^y in 
the exterior of the cylinders. 

Let Fx : En —> M be given by 

F,(77) := Yl ^l{v/N)Fxy{ij), 

yeS\{x} 

where each -q/N is thought of as a point in D and {9^ : y G S'\ {a;}} is the partition 
of unity established before. 

The following properties of F^ are helpful in the proof of Proposition 15.11 It is 
easy to check that 

FM = F^yil) for r;/7V G . (5.4) 
Indeed, if -q/N belongs to 3C^, Qy{r]/N) = 1 proving the identity claimed. On the 
other hand, li-q/N belongs to Df^, by definition of H, F^^{r]) — 1 for all z £ S\{x}, 
1 < J < — 1, so that Fxz{r]) = Fx{rj). By similar reasons, 

Fx = 1 on {rjeEN ■.rix>{l- 3e)iV} 

and Fx = Q on {r] £ En : Vx < 2eN} . (5.5) 

The minimum in definition (j5.3[) is introduced precisely to fulfill the second assertion 
in (|5.5|) . In particular, if ij/N e D^^ for some z G S then 

F4f]) = l{z^x}. (5.6) 
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Since H, as well as each Qy, is a smooth function, there exists a finite constant 
Ce, which depends on e through the definition of the smooth functions, but does 
not depend on > 1, such that 

m^\F,{a^^r^)-FM\ < § (5-7) 

for every z, w e S. 
Let 

3";^ := {iieEN ■■ilx + riy>N -In) , x^yeS . 
Clearly, 3^ — U^, x, y £ S and, for every N large enough, C Hf^ . Let 
3^ :— Uy^g\^^j3'^ . In what follows, the value of the constant may change from 
line to line, but will never depend on N. 

Lemma 5.2. For each x Cz S and every N > 1 large enough, 

Proof. By property (jS.Sp . we can restrict the sum in the left hand side to configu- 
rations 77 e En \ 3^ satisfying eN < rj^ < {I — e)N. So, by (|5.7p . the left hand side 
of the above inequality is bounded above by 

^ E ^^('?) ■ 

eN<ri^<(l-e)N 



This expression is bounded above by 



ZnsN^ ^ ^ ain) 

max{ ?7y : y } < — 2 — £jv 



which can be re-written as 



eN<t<il-e)N ^ ^ ^ ^ Ce£;iv-i.s\{x}(^Jv) 

By Lemma 13.21 for the expression inside braces, last expression is bounded above 

by 



By Lemma 13^21 once more and Proposition 1 2 . II we obtain the desired result. □ 

Fix a nonempty subset C. S and denote 5^ :— S \ ^ 0. We define the 
function Fgi : Ef4 ^ R as 

Fs^{ll) := E ^-(^) ■ 

xes^ 

Let us define the following subsets of i?jv 

Dl, := {7]eEN:Vx>N- Mn} , xeS, 
so that £^ C Dff, X ^ St,. It follows from (|5.5p that if 7/ G Df^ for some x ^ S then 
Fsi(ry) = l{xe5i} = ^^s^ (a^'"/?) , (5.8) 
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for every z,w € S and every N large enough. In particular, 

Fs^ e Cn{ U D%, y Dj,) . 

xes^ yes^ 

We shall use Fsi to get an upper bound for cap^( D^eS^ Df^jLly^s^^N) ■ 
We first claim that for any N large enough, 

Fsiicr'^^T]) = 1 = Fs^iv) foraUr/e |J 3""^^ and z,w £ S . (5.9) 

To prove this claim, Rx x ^ y in . By ([5?5|l . for 77/iV e ilf*', 

W = F^yir,) + Fy^r^) . (5.10) 

Recall from (|5.2p the enumeration of 5 defined according to the values of fxy. 
Let 01, . . . , Zk and wi, . . . ,Wi^ be such enumerations obtained from f^y and fyx, 
respectively. Since fxy + fyx — 1? we can choose the enumerations in such a way 
that Zn+i — Wi^-n, < n < k— 1. With this convention, an elementary computation 
shows that 

K-l 

F.yiv) + Fy^v) = ^{fxy(z,) - fxy(z,+i)}( F4(r7) + Fy--^{rj) ) . 
By (|5.ip . the previous expression is equal to 

K-l 

Claim (|5.9p follows from this identity and (|5.10p since 3^ C XL^^ for N sufficiently 
large. 

For each subset A C En and function F : En ^ M., let 

^ > > srlTj|(7(T7.lrl2;-7/;l-!^liT'""T7l — /^(ril!-'^ 



DNiF;A) E f^Niri)9ivMz,w){Fia^^v) ~ F{7^)} 

With this notation, Lemma [5^21 can be stated as 

DNiF.;EN\rN) < j;^^^:^^^^ yxeS. 
By Cauchy-Schwarz inequality, 

DNiFsi;EN\U,es^3%) < \S'\ Dn{F,; En \U,es^3lf) 

< \S^\Y, DN{F^;EN\rN) . 

Therefore, since £n T it follows from Lemma [5T^ that 

lim N"+'Dn{Fsi;En\U,^si3'n) = 0. (5.11) 

It remains to estimate DAr(i^5i ; U^ggi J|^). By definition of 3n, z G 5^, and by 
(EH), 

Dn{fs^; U a^) = i?Ar(i^si; U ^^') = E E ^A^(^si;3^') • 
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The last identity follows from (15. 8p and the relation 

^n""" n C IJ i?^ for all XI, e and 2/1,^26^^ 

Therefore, by ()5.5|) and ()5.4|) we finally conclude that 

U = E E (5-12) 

We now provide an estimate for each term in this sum. To derive this bound, in 
addition to the properties already imposed to the function we also require that 

sup{0'(u) : It e [0,1] } < l + (5.13) 

and 

sup { : u e [2e, 1] | < 1 + . (5.14) 

u — e 

The first requirement can easily be accomplished since (1 + ^/e) times the length of 
the interval [3e, 1 — Be] is strictly greater than 1 for e small enough. For (I5.14|) . it 
suffices that 0(u) < (u — e)(l + y^) for all u e [3e, 1] because (/> vanishes on [0, 3e]. 
Since (w — e)(l + ^/e) > 1 for u = 1 — 3e and every e small enough, it is possible to 
define a smooth function cj) satisfying (j5.14p without violating the other previously 
imposed properties. 

According to the above discussion, in what follows we suppose that e is an 
arbitrary number in (0, ep] for a suitably chosen cq > and that </> satisfies the 
additional properties (|5.13p and (|5.14p . 

Proposition 5.3. For any x,y E S , x y, 

limsupiV"+ii?jv(i^,,;a^y) < ilt^^^l!^ caps(x, y) y G ^4 • 

Proof. Let x = zi, Z2, ■ ■ ■ , z^^ = y be the enumeration established in the definition 
of Fxy, so that fxy(2;n) > fxy{zn+i), 1 < n < K — 1. Fix two different sites Zi ^ Zj 
in S with 1 < i < j < k. By definition of F^y, 

j-i 

n—i 

Thus, by the Cauchy-Schwarz inequality, the sum 

E f^N{v)9{Vzi)r{zi,Zj ){F^yia^-^'v)-F,y{v)Y (5.15) 

is bounded above by {{xy{zi) — {^y{zj)} times 

5](fxy(^„)-fxy(^„+l)) J2 ^^NirMv^M^^,^MP^yi^''''V) - F^yiv)}' . 

Performing the change of variables ^ = rj — d zi , the second sum above is less than 

m4z^)r{z,,z,)^ J2 ^R",(0-^^",(^ + f.J} , 
^^^s a(0 
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where := e En-i,s ■ ^x+^y > N — 2£n}- So far, we have shown that ()5.15|) 
is bounded above by mi,(zi)r(zi, Zj)iV"Z^^g{fxy (z;) — ixyizj)} times 

Fix some i < n < j. The second sum in the above expression may be re- written as 

E^^" ^ mj ^ m^,{xYm^,{y)^-'^-'^ 

^ a(C) ^ a(fc)a(iV - m - fc) 

m=0 CGB„,s\{^,a} £jV<fc<(l-2e) AT ^ ^ ^ ' 

{i/(dfe+i(C))-ff(4(C))}', (5.16) 

where 

4(0 ]^ + l{">2}^^, 6iV < fc < (1 - 2e)iV . 

To keep notation simple let (pk stand for (j}(dk{C))- By the Cauchy-Schwarz inequal- 
ity, the last expression is less than {m^,{x)mi,{y)Y^ N"^" times 

y y Hh^ y u'^{i-urdu / Z — du. 

Since m < 21^ then, for all N large enough, the last integral above is less than 



{0fe+i-(/.fe}( sup {0(7.)/(w-e)})'" < 
«e[o,i] 



The last inequality follows from assumptions (|5.13p and (|5.14p . Therefore, we 
conclude that (|5.16p is bounded above by 

{m,(x)m.(i/)}^^(l-fVi)^"+i mi 
Io.N^'^+^ ^ ^ a(C) ' 

which in turn is bounded by 

Zsna){m,{x)m.{v)Y''{l + ^ef''+^ 

^ -A- r a; r y /q, ^ 

Hence, we have shown that f|5.15p is bounded above by 

, 2 ( Z5r(«){i -f v^p^+i ^ ^^^^^ 

In a similar way we can get the same upper bound for (|5.15p if we suppose 
instead that j < i. The assertion of the proposition follows from this estimate and 
Proposition □ 



.(z.)r(z„z,-){f3.y(z.)-fxy(^,)r( Ir 7 rv T ){m.{^)^n.{y)r^ N 



We are now in a position to prove Proposition 15.11 Let Sl := C] Si,, 
S'^ n 5*^, and suppose they are both nonempty sets. Since 

c&pj^{£NiSl),eNiS^,)) < cap^( y D%, [j Dl) < i^jv(Fsi) , 
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it follows from (l5lT|) . (|5?T2l) and Proposition [531 that 

limsupAr"+icapjv(£Ar(5i),£Ar(52)) < caps(x,2/), 

after letting e J, 0. Theorem 12.21 follows from Proposition l5.ll and Proposition 14. II 

6. Proof of Theorem 12.41 

In [2], we reduced the proof of the metastability of reversible processes to the 
verification of three conditions, denoted by (HO), (HI) and (H2). The proof of 
condition (HI) is similar to the one presented in [3] for zero range processes on 
complete graphs. However, in the case where m is not uniform, some modifications 
are needed to handle sites not in 5'*. This is the only reason for which we have 
introduced the sequences b^iz), z £ S \ S^, and the respective condition in (|2.10p . 

The following notation will be used throughout this section. For each cc € 5*, 
let S En be the configuration with N particles at x and let £^ represent the 
set £iv(5'* \ {x}). 

Condition (H2) follows immediately from ()3.2p since ^n{£-n) ^1/'** for every 
X € 5**: 

MM = 0, yxeS.. (H2) 

Fix a configuration rj in £^. Since J2y^x % — and ry^ < fe^r, for z G 5 \ 5*, 
by the explicit form of /i^f, 

(6.1) 

> TT m.(z)''«(^) . 

Hereafter, Cq stands for a constant which does not depend on > 1 and whose 
value may change from line to line. To estimate the capacity, cap^({7y}, {^^}) we 
consider a path r]'^^^ , < j < p, from 77'-°-' = 77 to 77^^-* = ^f; obtained by moving to 
X, one by one, each particle. Since there are at most particles to move, we can 
take a path such that p < kIn- Let F be an arbitrary function in CAr({r7}, {^^}). 
By Cauchy-Schwarz inequality and the explicit expression of the Dirichlet form. 



,=0 7^0 A^Jv('7(-'')) 



Therefore, by (|6.1 



cap^({^},{c^}) > n "^4^)""^^^ 

The extra factor £n comes from the length of the path. Condition (HI) follows 
now from this estimate. Theorem 12.21 and condition (j2.10p : 

■ , T'^^^T;\^L:.^^^ = » ' ^ 5. . (HI) 
iv^oo mf^ggx {cap^({r/},{^^})} 
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Finally condition (HO) follows from Theorem 12. 21 as we show below. Denote by 
Rf^ (•, •) the jump rates of the trace process {jyf " ■ t > 0} defined in Section[2l For 

X ^ y in Si,, let 



rN{x,y) := ^ MwW^w('?,0 



By Lemma 6.8 in [2], 
IJNi£-N)rN{£.N,£.N) = i|cap^r(£^, £^) + Capjv(£Ar , £a,) 

- capjv(^£Ar({x,y}) , £jv(5'* \ {a;,?;})^ | . 
Therefore, by Theorem 12.21 since /iAr(£^) converges to k^-^ for all a; in 5*, 

hm iVi+"r^(£^,£^^) = Vx, y G 5. , x ^ y . (HO) 

This proves Theorem 12.41 as a consequence of Theorem 2.10 in [5]. 
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